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Abstract 

We theoretically study the electronic transport properties of Dirac fermions through one and 
double triangular barriers in graphene. Using the transfer matrix method, we determine the trans- 
mission, conductance and Fano factor. They are obtained to be various parameters dependent such 
as well width, barrier height and barrier width. Therefore, different discussions are given and com- 
parison with the previous significant works is done. In particular, it is shown that at Dirac point 
the Dirac fermions always own a minimum conductance associated with a maximum Fano factor 
and change their behaviors in an oscillatory way (irregularly periodical tunneling peaks) when the 
potential of applied voltage is increased. 
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1 Introduction 



Graphene is a single two-dimensional array of carbon atoms with a honeycomb lattice, which was 
discovered in 2004 [I]. This finding has been attracted an intensive attention from both experimental 
and theoretical aspects. In particular, the tunneling of Dirac fermions in graphene has already been 
verified experimentally p|, which in turn has spurred an extraordinary amount of interest in the 
investigation of the electronic transport properties in graphene based quantum wells, barriers, pn 
junctions, transistors, quantum dots, superlattices, etc. The electrostatic barriers in graphene can be 
generated in various ways [3|[4]. For example, it can be done by applying a gate voltage, cutting it into 
finite width nanoribbons and using doping or otherwise. Whereas magnetic barrier could, in principle, 
be realized with the creation of magnetic dots. As far graphene, results of the transmission coefficient 
and the tunneling conductance were already reported for the electrostatic barriers (#[9] and magnetic 



barriers 10-12 



The electronic band structure (energy dispersion relation) of graphene consists of two inequivalent 
pairs of cones with apices located at Brillouin-zone corners. The dispersion relation E = zizhvF\k\ is 



linear around the Dirac point {K, K') where vp — 9.84 x lO^m/s is Fermi velocity 13 . The presence 
of such Dirac-like quasiparticles is expected to induce some unusual electronic properties, which make 
difference with respect to two-dimensional electronic gas, such as the so-called Klein paradox [14] , 
anomalous integer quantum Hall effect 15-17] and observation of minimum conductivity 16 . The 
fact that in an ideal graphene sheet the carriers are massless, gives rise to Klein paradox, which 
allows particles to tunnel through any electrostatic potential barriers, that is the wavefunction has an 
oscillatory tail outside the electrostatic barrier region. Hence this property excludes the possibility to 
confine electrons using electrostatic gates, as in usual semiconductors. Thus to enable the fabrication 
of confined structures, such as quantum dots, we need to use another type of potential such as the 



scalar potential coupling 18 



Theoretical investigations have been widely performed to clarify the resonant-tunneling features 
using mostly barriers of the rectangular forms. The reasons because the corresponding models are 
so simple to have an advantage for numerical calculations. However few works studied tunneling 



effect with barriers of the potential slopes as a result of externally applied field 19-21 . One of 



them is the trapezoidal double barrier structure, which was investigated to study the effect of the 
potential disturbance at the interfaces of the graphene cheet [23] . In the same spirit, we consider 
another problem based on single and double triangular barrier structures. Our model is possibly 
applied to the resonant tunneling diodes of which the barriers are formed by delta doping in the 
future and is also a step outward from a rectangular form from the other point of view. We ensure the 
confinement of Dirac fermions in the y-direction by using infinite mass confinement, which requires 
infinite mass at the boundary of the y-strip and results in a specific quantization of the y-component 
of the momentum [l8] . The effects of the well width, barrier height and barrier width on the transport 
properties are systematically studied through numerical calculations. As long as the applied potential 
is increased, the number of the minimum conductance associated with maximum Fano factor increases 
as well. This result makes difference with respect to that of rectangular barrier where there is only 
one minimum and one maximum [9]. We conclude that it is relatively more easily for Dirac fermions 
to tunnel through a triangular barrier in a graphene sheet rather than rectangular one. 
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The outline of the paper is the following. In section 2, we set our theoretical model by giving the 
appropriate equation describing Dirac fermions in graphene and choosing the convenient configuration 
for the triangle double barrier structures as depicted in Figure [ij In section 3, we expose the exact 
analytical solution to solve the Dirac equation in each regions of the structures, which resulted in giving 
the corresponding eigenvalues and eigenspinors. Tunneling probabilities are calculated in section 4 as 
a functions of different parameters such as the fermion energy, static electric field and incident angle. 
These are done by matching spinors in different interfaces and using the transfer matrix techniques. In 
section 5, we discuss the transport results corresponding to single and double barriers separately. The 
obtained results show characteristic oscillations associated with tunneling resonances as a function of 
the fermion energy and the static electric field. We conclude our work in the final section. 



2 Theoretical formulation 

We consider a system of massless Dirac fermions through a strip of graphene characterized by the 
length Lb and width w in the presence of a double triangular barriers. In the systems made of 
graphene, the two Fermi points, each with a two-fold band degeneracy, can be described by a low- 
energy continuum approximation with a four-component envelope wavefunction whose components are 
labeled by a Fermi-point pseudospin = ±1 and a sublattice forming an honeycomb. Being a zero-gap 
semiconductor, the quasiparticle motion can be described by the massless Dirac like equation 

[vF^-p + V{x)]Tp{x,y) = e'tp{x,y) (1) 

where vp — 9.84 x W^m/s is the Fermi velocity, p = —iK^ is the momentum operator (hereafter 
vp = h = \), a = (fja;, (Ty) are the Pauli matrices, e = hvF\k\ being the energy of the incident particle. 
A triangular double barrier configuration is depicted in Figure [l] with all parameters, which requires 
two kinds of width: the right and left sides of the barrier. Therefore the dependence of the various 
parameters can be considered as shown in the potential V{x) configuration 



V{x) 



where we have set Fi = F2 = F^ = and F^ = are the strength of the static electric 
field in each regions. 

Our system is supposed to have finite width w with infinite mass boundary conditions on the 
wavefunction at the boundaries y = and y = w along the y-direction [6l[7l [Io[jT2|[l8] . This boundary 
conditions result in a quantization of the transverse momentum along the y-direction, which is 

ky = kn = ^(^n+^^ , n = 0,l,2---. (3) 

One can therefore assume a spinor solution of the following form tljj{x,y) = (0j^(a;), (/>j_ (x))*e*'^J'^ 
where j = 1 for x < a, 2 for x £ [a, 6], 3 for x £ [6, c], 4 for x £ [c, d], 5 for x £ [d, e], 6 for x £ [e, /] and 
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Figure 1: The parameters of a triangle double barrier structure. 



7 for X > f denotes the different space regions. Thus our problem reduces to an effective ID problem 
whose Dirac equation can be written as 




dx + 




0. 



(4) 



-i(t)j_{x) 

An electron which impinges from x = —oo on the quantum barrier is partially reflected, partially 
transmitted at the interface x = a. Inside the barrier regions [a, c] and [c?, /], the eigenstates is a 
combination of the parabolic cylinder functions waves. For x > f, the carrier is also partly transmitted 
and escapes towards x = oo with a wavevector +kx. The electric potentials Fi, F2, F3 and F4, being 
uniform along the y-direction, the y-component of momentum is conserved throughout the regions. 
Due to the space dependence of the potential V{x) we make the following transformation on our spinor 
components to enable us to obtain Schrodinger like equations for each component, Xj+ = 5 {4'j+ + 4>j-) 
and Xj- = ^ {4'j+ ~ 4>j-)^ which obey the coupled stationary equations. These are 

d 



dx 



Xj± (x) ± i {V{x) - e) Xj± (x) =F ikyXj^ 



0. 



(5) 



Each spinor component Xj± can be shown to satisfy the following uncoupled second order differential 
equation 



Xj± (x) 



±i^V(x) + [V(x) 
dx 



Xj±{x) = 0. 



(6) 



At this stage, we point out that our effective 2D massless Dirac equation Q is equivalent to a massive 
one with an effective mass m equal to the transverse quantized wave vector ky, i.e. m = ky = kn- For 
this purpose, we consider a unitary transformation, which enable us to map the effective 2D massless 
equation into a 11? massive Dirac equation. Such a unitary transformation does not affect the energy 
spectrum or the physics of the problem. We choose a rotation by 7r/4 about the y-axis, U = e^i"^ 
and thus the transformed Hamiltonian and wavefunction read 



V{x) — £ + ky 



A. 

dx 



dx 



V{x) 




^j+,-{x) = U^j. ix). 



(7) 
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which is identical to a ID massive Dirac equation with an effective mass m* = ky. This shows clearly 
how to derive the dynamical mass generation via space compactification [22] from our model 

3 Exact solution 

After solving the differential equation Q, It turns out its solution in regions x < a, [c,d] and x > f 
are given by 

Mx) = ( ^ ]e"^^- + rj _] )e-'=^- (8) 

A.( \ { ^4 { ^^(^ ^ 

04(x) = a„4 , +Pn4 _, , 

V % (^) / V ^4 / 

= tn\ ^ )e*'-^ (9) 

where r„ and i„ are the reflection and transmission amplitudes, respectively, n is labeling the modes, 
the functions (x) and v^{x) are u[{x) = (x) = u^{x)/z = —v^ (x) = e*'^"'^. The wavevector 
kx = ^ e"^ — k^ and the complex number Zn^k^ is defined as 

e) — , (10) 

where the transversal momenta is quantized as shown in ([s]). Note that this quantization is the 
result of the infinite mass boundary conditions mentioned previously on the wavefunction along the 
y-direction. In the case of \ky\ > \e\, the waves are evanescent (bound states) outside and inside the 
quantum barrier and thus the imaginary wavevectors associated with the evanescent waves are given 



by kx = i^ky — e"^. Since we are interested by the transmission of relativistic particles (continuum 
scattering states), thus we disregard the bound states which correspond to imaginary kx • 

The solution of ^ in the quantum barrier (region [a,b]) can be expressed in terms of the parabolic 
cylinder function Dy[x) as 

X2+[x) = a„2-D!.„i_i {qi) + Pn2D^u„i {-ql) (11) 

where Vni = '^^"2Vi^^ , qi = \J (a-b)Vi ^^'^^'^^'^^ ~^ -^i)' ~ —aVi + (a — h)e, the parameters a„2 and 
/3„2 are constants. Now substituting ( [IT| ) into ^ to get the second component of X2{x) 

X2.(x) = an2 , ^2{a-b)Vie-'^/^D,^,{q,) (12) 

(a - b)kn 

+ V^2(^^6yi4e*"/^Z)i-.,, {-ql) 2{-Vix - Ei)D.,„, {-ql)]. 

The components of the spinor solution of the Dirac equation ^ in the region [o, b] can be obtained 



from (11) and (12) with (f)2^{x) = X2+ + *X2_ and (f)2-{x) = X2+ — iX2-- These give 



fcw="n.ri!i + A. (13) 

V U2 [X) I \ V2 [x) I 
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where the functions uf{x) and vf{x) read as 



±1 



vf{x) 



a — b k„ 
1 



(a - 6)3/2fc„ ' 

+ Va-6 {b{±2is - kn) - a{±kn + 2i{Vi - e)) T ^iVix) D^y^^ {-ql)] 



(14) 



Similarly, the solution of ([g]) in the region [6, c] takes the form 

X-i+{x) = an3-Di.„2 (92) + l3n3D-y„2-i (-12) 



(15) 



, ikl{h-c) 

where Vn2 = 2V1 ' ^2 
X3(a:) is given by 



^e^'^^^{Vix + £"2), £'2 = — cVi + (c — b)e. The other component of 



«n2 ,, \, [-V2(5-c)yie-^"/^D,,+i (gs) + 2(^1^ + ii;2)i?.„, (92)] 

[0 — CjKn 



+ /3n2 ,, \, V2(6-c)yie--/^Ij-.„,(-g2*). 



(16) 



Combining (15) and (16) in similar way to (f>2{x), we obtain the eigenspinor solution of the Dirac 

v+(x) \ 



equation Q in the region [b, c] 



u+(x) 



(17) 



where we have set 



uf{x) 



1 



-[Tv^(6-c)e-/4Z),„,+i((72) 



(6 - c)3/2A:„ 

+ (6(T2ie + kn) - c{kn + 2i{±Vi ^ e)) ± 2iVix) D^^^ (92)] 

±A/^^e— /^Z)_.„, {-ql) + Z)_,„,_i {-ql) 



(18) 



V b — ckn 

Finally, note that the general solution of equation ^ in regions [d, e] and [e, f] can be obtained by 



interchanging a^d, b^e, c^f and Vi — )• V2 in the equations (14) and (18). The coefficients r„, 
tn, Cinj and (inj {j = 2, 3, 5, 6) can be determined by matching wavefunction at different interfaces. 

4 Transport properties 

The transmission coefficient is determined by imposing the continuity of the wavefunction at the 
interfaces between regions. This procedure is most conveniently expressed in the transfer matrix 
formalism. Here we directly use this approach and refer the reader to references [24y25| for a detailed 
discussion. The transfer matrix defined by 



M 



Mil M12 
M21 M21 



(19) 
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relates the wavefunction on the left side of the barrier structure {4>x<a) to that on the right side ((pxyf)- 
We can then construct 2x2 matrices in each region, whose columns are given by the spinor solutions, 
such as for regions (1,4,7) 



wi{x) = W4{x) = wrix) = * -ik X ] 



and for j = 2, 3, 5, 6 

w,{x)=( ""^^'^ ^^"0 (21) 

These matrices play the role of partial transfer matrices and allow to express the continuity condition 
of the wavefunction at each interface. Note that the components of matrices W5{x) and wq{x) can be 
obtained by interchanging a ^ d, b ^ e and c — )• / in matrices ^2(2;) and ^3(2;), respectively. After 
straightforward algebra, we get the transfer matrix as function of different boundaries 

M = w^\a)w2{a)w^\b)w3{b)w^\c)wi{c)w^\d)w5{d)w^He)we{e)w^\f^^^ (22) 

and the relation which expresses the continuity of the wavefunction is then given by 



Solving equation (23) for the transmission amplitude t„ of the nth wave mode through the barrier, we 
get the transmission probability T„ as 

2 1 

^M^nF' ^^^^ 

Based on [26| , we give a review about the shot noise. Indeed, the conductance of a single trans- 
mission channel can be written as ^ 

G = g^T (25) 

where g is the degeneracy (spin and valley) of the system and T the electron transmission probability. 
When the system is biased, shot noise appears due to discreteness of charge 27 and these current 
fluctuations for a single channel are given by 

{{6lf)=2e{I){l-T). (26) 

The total noise power spectrum for a multichannel conductor is then obtained by summing over all 
transmission eigenchannels: 

5/ = ^'^£~V„(l-r„). (27) 

n=0 

In the limit of low transparency T„ <C 1, 

Si = 5poisson = ^" = (28) 

n=0 
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defining a Poissonian noise induced by independent and random electrons like in tunnel junctions [27] . 
The regular way to quantify shot noise is to use the Fano factor F which is the ratio between the 
measured shot noise and the Poissonian noise: 

5poisson 2e(I) j:^^-^^t^ ■ 

Then, for a Poissonian process = 1 at small transparency (T„ — )• 0), while F = in the ballistic 
regime (i.e. when T„ — )• 1) and F = 1/3 in the case of a diffusive system [9,28-30 



In graphene, it has been theoretically concluded that transport at the Dirac point occurs via 
electronic evanescent waves [9,31 . Tworzydlo et al. [9] used heavily-doped graphene leads and the 



wavefunction matching method to directly solve the Dirac equation in perfect graphene with length 
Lb and width w. They found that for armchair edges, the quantization condition of the transverse 
wave vector is defined by 

(n + a) , ^ 

ky,n = ^ 30 

where a = or 1/3 for metallic and semiconducting armchair edges, respectively. At the Dirac point, 
the transmission coefficients are given by j9j 

Tn = 7 (31) 

cosh I TT (n + a) - 



Consequently, graphene has a similar bimodal distribution of transmission eigenvalues at the Dirac 



point as there is in diffusive systems 32 ,33 . In the limit of w/Lb — )• oo, the mode spacing becoming 
small and one can replace the sum over the channels by an integral over the transverse wave vector 
component ky to obtain the conductivity and the Fano factor for a sheet with metallic armchair edge [9] 

Lb Ae^LB r dky Ae' 
w h w Jq cosh [ky Lb) ^rft 

77 _ Sn=0 '^n{^ — Tn) 1 , , 

''^Dirac - ^N-lrr " q' ^"^"^^ 

Z^n=0 ^ri -3 

In summary, we will study the above quantities for the present system in terms of our findings 
and compare with already published works. In fact, the conductivity and the Fano factor of Dirac 
fermions through one and double triangular barriers in graphene will have a variate and different from 



with respect to the results presented in (32) and (33). 



5 Results and discussions 

For a better understanding of the obtained results so far, we numerically analysis different physical 
quantities in terms of the system parameters. To underline their behaviors, we trait single and double 
triangular barriers, separately. 

5.1 Single barrier 

We start our discussion by studying the transmission probability and shot noise for the Dirac fermions 
scattered by a single triangular barrier potential. We implement our previous analytical approaches to 
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a graphene system subject to a single triangular barrier potential of strength Vi and V2 = 0. We will 
see that the transmission coefficient has a rich information about the electronic transport properties 
of the Dirac fermions through a triangular barrier structure. 

The variations of the calculated transmission coefficient T in terms of the incident electron energy 
e and applied voltage Vq is displayed in Figure [2] for different values of the barrier widths Lbi, barrier 
widths right side Lir, barrier heights Vq and incident electron energy e. From Figure [2j one can 
see that for certain values of Lbi (Figure [2^), Vq (Figure [2j3) and Lir (Figure [2|:) , the transmission 
resonances appear in the triangular potential for m* < e < Vq + 2m* and vanishes for both conditions 
{e > Vq, Lir ^ Lii). We note that the intensity of resonances increase as long as Vq and Lbi increase, 
which allow for emergence of peaks in the T shape. One can see that T always starts from the energy 
corresponding to ky = m* , with m* is the effective mass of the ID Dirac fermion. The zone when 
we have the energy such as e < m* corresponds to the forbidden zone. It is important to note that 
the resonant energy depends strongly on the barrier height and width. For L^i = 4 we have the 
transmission resonances independtly of the value taken by the applied potentail Vq as long as Vq < 
While for Vq > e — 2m*, the resonances decrease sharply until reach a relative minimum and then 
begin to increase in an oscillatory manner. 




Figure 2: Transmission coefficient T for the Dirac fermion scattered by a single triangular barrier potential 
with Vi = Vq and V2 = 0. (a), (b) and (c) as a function of incident electron energy, (d) as a function of 
applied voltage. 

Figure |3] is showing the transmission coefficient T as function of the electron incident angle 4> for 
e = 2 and different values of (Vq, Lbi, Lir). We see that the perfect transmission occurs at different 
angles and vice versa. It is observed that, the transmission is always total for a normal incidence angle. 
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For Vo = 5 one can observe that T is not zero for some values of the barrier width. In particular 
it shows up two peaks at incident angles (f) = ibg-g|^ and (p = itg^^gg for each value of the barrier 
widths L^i = 10 and Lbi =4, respectively. The transmission resonances always appear for the case 
of the barrier width only of the right side Li^. is equal to the barrier width only of the left side Ln, 
i.e. Lir = Lii, while disappear otherwise {Lir ^ Ln). 




(c) 


m 











Figure 3: Transmission coefficient T as a function of electron incident angle (p for the Dirac fermion 
scattered by a single triangular barrier potential with Vi = Vq, V2 = and e = 2. 



In what follows we discuss the conductivity a and Fano factor F behaviors to underline what 
makes difference with j9 , 16 ,34 . Indeed in Figure [4| 




Figure 4: The electrostatic potential dependence of the Fano factor and the conductivity for the Dirac 
fermion scattered by a single triangular barrier potential with Vi = Vq, V2 = 0, and e = 2. 
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we plot a (in units of Ae^/irh) and F in terms of the electrostatic potential Vq for Lbi = 2 and e = 2. 
It is interesting to note that a corresponding to our system is showing some differences with respect 
to that for an ideal strip of graphene [o] , which supports perfect transmission regardless of the barrier 
height (Klein tunneling j|3|). It is obvious to observe that the conductivity and Fano factor change 
their behavior in an oscillatory way (irregularly periodical tunneling peaks) when we augment the 
potential of applied voltage. One can see that as long as Vq increases, the number of minimum of a 
increases as well but the associated number of maximum of F decreases. This effect is different from 
that obtained in |9| where there is only one minimum conductance 4e^/(7r/i) at the Dirac point and 
for a geometric factor w/Lb = 5, which corresponds to one maximum Fano factor 1/3. In contrast, in 
our case for the same factor (w/Lb = 5) the minimum conductivity 0.387 appears around Vq = 3.661 
where the associated maximum Fano factor is 0.315. More importantly, for two values of Vq like 3.661 
and 7.616 one can see that the minimum conductivity increases from 0.387 to 0.479. Therefore, we 
conclude that it is relatively more easily for the Dirac fermions to tunnel through a triangular barrier 
in a graphene sheet rather than rectangular one. It should be pointed out that the nonzero minimum 
conductance, as shown in Figure |4| may due to the conservation of pseudospin and the chiral nature 
of the relativistic particles in the graphene nanoribbon. 



5.2 Double barriers 



In this section we implement our previous analytical approach to a graphene system subject to a 
double triangular barrier potentials Vi and V2 so that the resulting static electric field strengths are 

Note that there are various parameters involved such 



^1 = 1^, 



^, F3 = ^ and F. 



.Y2 



as well width, barrier height and barrier width, these will offer different discussions about transport 
properties in the present configuration of potential. In particular, Figures [5] and [6] show the transmis- 
sion coefficient in terms of the incident electron energy e and applied voltage Vi = V2 = Vq for both 
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Figure 5: Transmission coefficient T for the Dirac fermion scattered by a double triangular barriers potential 
with Vi = V2 = Vq, Lw = and Lbi = Lb2- (a), (b) and (c) as a function of incident electron energy, 
(d), (e) and (f) as a function of applied voltage. 
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cases Liij = and 7^ 0, with is the interbarrier separation (weh width). In Figure [sj one can see 
that contrary to single barrier (e.g. Figure 2c) the transmission resonances always appear for double 
triangular barrier case. Clearly, the intensity and width of resonances as well as the condition for the 
existence of resonances depend on the static electric field strengths {Fi, F2, F3, F4) and barrier widths 
{Lbi, Lb2)- The intensity of resonances increases and decreases as long as the strengths IF2I, 
I-F3I, IF4I) decrease and widths {Lbi, Lb2) increase, respectively. 
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Figure 6: Transmission coefficient T as a function of incident electron energy for the Dirac fermion 
scattered by a double triangular barriers potential with Vi = IV2I and ky = 1. (a) the width of the well 
was varied, (b) the barrier height is varied, (c) the barrier width only of the outer sides is varied, (d) the 
height of the barrier and the well is varied. 



In Figure [6] we consider the same behavior of transmission as before but with 7^ 0. Compared to 
Figure [5] for Vi = V2 = W, one can conclude that the intensity of resonances depends strongly on Lw 
In addition for Vi = V2 there are several peaks showing transmission resonances those correspond to 
the bound states and no resonances exist for Vi = —V2. Figure [6^ presents three potential profiles 
and their transmission coefficients in terms of e for Lw = 1,2,3, Lbi = Lb2 = 2 and Vi = V2 = 10. 
The results show that as long as the well width Lyj increases the resonance peak structures become 
sharpened. We consider in Figure [6|3 the case of three values of barrier height (Vi = V2 = 5, 10,20) 
for Lbi = Lb2 = and L^ = 3. As the barrier height is enhanced, the transmission resonance shifts 
and the width of the resonances increases. The case of the double barriers is more important when 



11 



the outer slopes of the barriers are varied with a fixed barrier height as shown in Figure [G]:. We notice 
that the weh region and the inner slope region are not changed when the parameters {Ln = L2r) vary. 
As long as the barrier becomes thick and high the resonance shifts toward the higher and the peak 
is more sharpened. Similarly to one triangular barrier, no resonance exists for the incident electron 
energy higher than Vi + 2m* and the zone e < m* is a forbidden zone. Except that when Vi = —V2 
no resonances exist for the following cases (1 < e < 9, Vi = —V2 = 5), (1 < e < 5, 14 = — V2 = 10) and 
(1 < e < 4.5, Vi = —V2 = 20), which are clearly shown in Figure [6jl. 

We represent in Figure [7] the transmission coefficient versus the incident angle with the same 
parameters as in Figure |3] for the Dirac fermion scattered by double triangular barriers potential with 
the interbarrier separation = 0. By contrast with the case for the Dirac fermion scattered by a 
single triangular barrier potential we conclude that the transmission resonances still always exist. The 
comparison between these two types of potential shows that for double barriers of strength Vi = V2 = 5 
we have three peaks with two peaks at incident angles 



±- 



TT 



vr 



"4.8021 
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and one peak at 



±7r 



(34) 



(35) 



for each value of the barrier widths L^i = 10 and L^i = 4, respectively. Even though the barrier 
width only of the right side Lir is different to the barrier width only of the left side Ln for double 
barriers, the transmission resonances always appear contrary to the Dirac fermion scattered by a single 
barrier. We observe that decreasing the barrier width only of the right side the transmission coefficient 
takes relevant values for a wider set of incident angles. 



(a) 




1 = 5 




Figure 7: Transmission coefficient T{4)) for the Dirac fermion scattered by a double triangular barriers 
potential with Vi = V2 = Vq, = and e = 2. 



As far as the conductivity and Fano factor behaviors for double barriers are concerned, we notice 
that the shot noise is characterized by the maximum of peaks at the minimums of conductivity and 
minimum of valleys at the maximums of conductivity. The role of the interbarrier separation L^, 
resulted in increasing peaks of shot noise and lowering the current of valleys as shown in Figure [8| 
One can see that the value F = 1/3 for the Dirac fermion scattered by double barriers is reproduced in 
the case where the barrier widths L^i = = 2, the interbarrier separation L^^, = and the applied 
voltage Vq near the 2e. 
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Figure 8: The electrostatic potential dependence of the Fano factor and the conductivity for the Dirac 
fernnion scattered by a double triangular barriers potential with Vi = V2 = Vq, and e = 2. 



6 Conclusion 

We have analyzed the behavior of Dirac fermions in graphene submitted to electrostatic potential 
of triangular type. By solving the eigenvalue equation we have obtained the solutions of the energy 
spectrum in terms of different physical parameters involved in the Hamiltonian system. Using the 
continuity of the wavefunctions at the interfaces between regions inside and outside the barriers, we 
have studied the transport properties of the present system. More precisely, using the transfer matrix 
method, we have analyzed the corresponding transmission coefficient, conductivity and Fano factor 
for single and double triangular barriers. 

It has been shown that the Dirac fermions scattered by single triangular and double triangular 
barriers own a minimum conductivity associated with a maximum Fano factor. We have noticed that 
the Dirac fermions can tunnel more easily through a barrier in the triangular forms rather than in the 
rectangular one. On the other hand, the behavior of the conductivity and the Fano factor in terms of 
the applied voltage showed irregular periodical oscillating for triangular potential. 

We have noticed that the resonant energy is influenced by the barrier width. Indeed, when the 
barrier becomes thick and/or high, the resonant peak becomes sharpened and shifted to the higher 
energy. Even if the thickness and the height of the barriers are constant, the form of the triangular 
barrier affects the resonant energy. On the contrary, the triangular double barriers structures is less 
sensitive to the well width compared with a rectangular double barrier structure. Therefore, we have 
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concluded that it is relatively more easily for the Dirac fermions to tunnel through a triangular barrier 
in a graphene sheet rather than rectangular one. These results may be helpful to deeply understand 
the transport in the nanoribbons and design the graphene-based nanodevices. 

We close by mentioning that the obtained results can be extended to deal with other issues related 
to graphene systems. Indeed, one may think to study the transport properties of Dirac fermions 
scattered by periodical potentials and other types. Another interesting problem, what about a gener- 
alization of the obtained results to bilayer graphene and related matter. 
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